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We consider spin-boson models composed by a single bosonic mode and an ensemble of N 
identical two-level atoms. The situation where the coupling between the bosonic mode and the 
atoms generates real and virtual processes is studied, where the whole system is in thermal equi- 
librium with a reservoir at temperature /3 . Phase transitions from ordinary fluorescence to 
superradiant phase in three different models is investigated. First a model where the coupling 
between the bosonic mode and the j — th atom is via the pseudo-spin operator a/U is studied. Sec- 
ond, we investigate the generalized Dicke model, introducing different coupling constants between 
the single mode bosonic field and the environment, g\ and #2 for rotating and counter- rotating 
terms, respectively. Finally it is considered a modified version of the generalized Dicke model with 
intensity-dependent coupling in the rotating terms. In the first model the zero mode contributes 
to render the canonical entropy a negative quantity for low temperatures. The last two models 
presents phase transitions, even when only Hamiltonian terms which generates virtual processes 
are considered. 
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1 Introduction 

The implementation of new techniques and ideas has lead to increase the interest in spontaneous 
radiation and collective effects of spin-boson models in free space and cavities [1] [2] [3] [4] . This 
has been a field of very attractive research, where experimental and theoretical progress have 
emerged, which can be useful in implementing concepts of quantum information and quantum 
computation in real physical systems [5]. The purpose of the present paper is to investigate 
different spin-boson systems, and the physical consequences that follow from assuming first that 
in interaction Hamiltonian of the models we are also considering the counter-rotating terms. In 
the interaction picture these terms describe non-resonant processes in which the two-level system 
and the bosonic field are excited or de-excited simultaneously. Second we study a interaction 
Hamiltonian that generates quite particular virtual processes. We consider in different models the 
possibility that virtual processes generate a quantum phase transition and also a phase transition 
at finite temperature from ordinary fluorescence to a superradiant phase characterized by the 
presence of a condensate. 

We consider three different models, in which we assume that a single mode of a bosonic 
field interacts with an ensemble of N identical two-level atoms. The whole system is in thermal 
equilibrium with a reservoir at temperature f3~ l . First, we study a modified version of the model 
discussed by Chang and Chakravarty, Legget and others [6] [7] [8], which has been used to analyze 
dissipation in quantum computers [9] [10]. Owing to the coupling between the two- level systems 
and a bosonic mode, and assuming that its intensity is generic, the zero mode contributes to 
render the canonical entropy a negative quantity for low temperatures. After dealing with this 
problem, the partition function is seen to be analytic for all temperature, and therefore there is 
no phase transition in the model. Second, we study the generalized Dicke model [11] [12] [13] 
[14] . We introduce different coupling constants between the single mode of the bosonic field and 
the ensemble of iV atoms, g\ and g2, for rotating and counter-rotating terms, respectively. In a 
situation where only virtual processes contribute, the generalized Dicke model present a second 
order phase transition from the ordinary fluorescent to the superradiant phase respectively, at 
some critical temperature f3~ l [15] and also a superradiant phase transition at zero temperature 
[16] [17], i.e., a quantum phase transition. We finally study the generalized Dicke model where the 
coupling between the bosonic mode and iV two-level atoms is intensity-dependent [18] [19] [20]. 
The intensity-dependent contribution appears only in the real processes. In the last model we are 
showing the presence of the quantum phase transition. 

The physics of superradiance is well known [21] [22] [23] [24] [25]. Let us suppose N identical 
two-level atoms prepared in the excited state. For dilute atomic systems, where there is no coupling 
between the atoms, the N atoms radiate incoherently, where the radiation rate is proportional on 
the number of the atoms N. Since the non-decay probability is exponential in time, the intensity 
of emitted light has the time dependence I(r) = Iq e ^~ ', where r e is the characteristic time for 
the spontaneous emission. Other characteristic of the ordinary fluorescence is that the radiation 
pattern of the atoms is isotropic. For an ensemble of atoms in a volume small compared to 
the emission wavelength, they start to radiate spontaneously much faster and stronger than the 



ordinary fluorescence case, where the radiation rate becomes quadratic-dependent on the number 
iV of atoms. Other important characteristic of this cooperative process is that the emission has a 
well defined direction depending upon the geometry of the sample. 

We are using path integral approach with the functional integration method to investigate the 
thermodynamic of the models, which is given by the analytic properties of the partition function 
in the complex /3 plane [26] [27] [28]. To study the nonanalytic behavior of the partition function 
of the models using functional methods two steps are mandatory. First, it is necessary to change 
the atomic pseudo-spin operators of the models by a linear combination of Grassmann Fermi fields 
to define fermionic models. Second, the thermodynamic limit (N — > oo), where N is the number 
of two-level atoms must be taken . For the first two models that we are interesting, the coupling 
between the pseudo-spin operators used to describe the N two-level atoms and the bosonic mode 
is linear. Consequently the path integral describing the ensemble of N atoms is Gaussian and the 
integration over the degrees of freedom describing the atoms can be performed exactly. 

A summary of results obtained in the present literature in the models is in order. The advantage 
of the first model is that it allows an exact analytic solution, and presents destruction of quantum 
coherence without decay of population. Concerning the thermodynamics of the second spin-boson 
model, i.e., the Dicke model, an important result was obtained by Hepp and Lieb [29]. These 
authors presented the free energy of the model in the thermodynamic limit. For a sufficiently large 
value for the coupling constant between the N two-level atoms and the single quantized mode of 
the bosonic field, there is a second order phase transition from the normal to the superradiant 
phase. Later, without assuming the rotating-wave approximation and by using a coherent state 
representation, the study of the stability of the model with an infinite number of bosonic modes 
was presented [30]. The study of the phase transitions in the Dicke model was presented also 
by Wang and Hioe [31], where some of the results obtained by Hepp and Lieb were reobtained. 
The generalized Dicke model, where the counter-rotating terms are also present in the interaction 
Hamiltonian, was investigated also by Hioe [12], Carmichael et al. [13], Duncan [14] and Li et al. 
[32] . Hioe studied the thermodynamics of the generalized Dicke model with two different coupling 
constants using also the coherent states. Li et al. pointed out that the A 2 term has been neglected 
in many papers, since the presence of such term leads to non occurrence of the quantum phase 
transition, as was proved by Rzazewski and Wodkiewicz [33] [34]. Carmichael et al. claim that 
in the limit N — »■ oo, the thermodynamical properties of the model are obtained by simply using 
the expression obtained from the model with only the rotating terms and doubling the coupling 
constant. As we will see, the thermodynamical properties of the generalized Dicke model deserve 
a more careful analysis [15]. 

More recently a bosonization procedure was employed to study the phase transitions in the 
generalized Dicke model. Employing a Holstein-Primakoff mapping [35] [36] [37], Emary and 
Brandes [38] [39] were able to express the generalized Dicke model in terms of a two mode-bosonic 
field. These authors discussed the relation between the quantum phase transition and the chaotic 
behavior that appear in the model for finite N, where the energy level-spacing statistics changes 
from Poissonian to one described by the Gaussian ensemble of the random-matrix theory [40] [41]. 
This chaotic behavior was discussed also in the Jaynes-Cummings model [42] by Graham and 



Hoherbach [43] [44] and Lewenkopf et al. [45] , in the situation where the counter- rotating terms 
are present in the interaction Hamiltonian [46] [47], since the seminal paper of Milonni et al. [48]. 

At this point, some remarks concerning the practical realization of the generalized Dicke model 
in the laboratory is in order. Dimer et al. [49] pointed out that it remains a challenge to provide a 
physical system where the counter-rotating terms are dominant. In the Jaynes-Cummings model, 
Cirac et al. [50] discussed the possibility of controlling the relative importance of the counter- 
rotating terms using a ion trap. The question that arises now is if there is some model where at 
some range of the physical parameters, the counter- rotating terms are dominant. It is possible to 
show that in the generalized Dicke model with intensity-dependent coupling in the rotating terms 
the relative importance of the contributions from the rotating and counter-rotating terms can be 
controlled by changing the temperature of the thermal bath. Concerning the intensity-dependent 
coupling model, Buck and Sukumar [18] [19] showed that the Heisenberg equations of motion can 
be solved exactly and the behavior at finite temperature was analyzed. Finally, Buzek [20] using 
the rotating-wave approximation also studied the same model, by presenting the time evolution 
of the model. 

This paper is organized as follows. In section II the path integral with functional integral 
method is applied to study the thermodynamics of a modified version of the model discussed by 
Chang and Chakravarty, Legget and others. In the section III we study the thermodynamics of 
the generalized fermionic Dicke model. In section IV we repeat our analysis for a generalized Dicke 
model by including the intensity-dependent coupling in the terms describing resonant processes. 
Conclusions are given in section V. In this paper we use the terms environment and reservoir for 
a system with a finite number of degrees of freedom and a system with an infinite denumerable or 
not degrees of freedom, respectively. In the paper we use ks — c — % — 1. 



2 The functional integral for the modified fermionic Chang 
and Chakravarty model 

Let us consider a bosonic system S, with Hilbert space 7i^ which is coupled with an ensemble 
of N two- level atoms, with Hilbert space HS B \ Let us assume that the complete system is in 
thermal equilibrium with a reservoir at temperature f3~ l . Let us denote by H$, Hb and Hi are 
the Hamiltonians of the bosonic field, the free N two-level atoms, and the interaction between 
both systems, respectively. The Hamiltonian for the total system can be written as 

H = H S ® Ib + Is ®H B + H T , (1) 

where Is and I B denotes the identities in the Hilbert spaces of the bosonic field and the ensemble 
of N atoms. 

The aim of this section is to study the analytic behavior of thermodynamical quantities, i.e., 
whether or not the system exhibits a phase transition from normal to superradiant phase at some 



critical temperature characterized by the presence of a condensate in a model similar to the one 
introduced by Chang and Chakravarty, Legget and others. Chang and others describes a system 
of one two-level atom coupled to a reservoir of harmonic oscillators, where the coupling between 
the reservoir and the atom is done via the pseudo-spin operator a z . They studied finite-time 
radiative processes in the model by assuming that the environment is in thermal equilibrium. A 
generalization of this model can be achieved introducing N identical two-level atoms and also an 
arbitrary mode-dependent coupling constant. Since [H,aZ\] = 0, this model allows for an exact 
analytic solution. The Hamiltonian of this generalized model reads 

Q N I N _ 

z 3=1 k V iv j=1 k 

where g^ describes the coupling between the j-th two-level atom with the reservoir. As usual, 
we are using the pseudo-spin operators at-., a7-s and of^ which satisfy the standard angular 
momentum commutation relations corresponding to spin | operators. We are also shifting the 
zero of energy \{oj\ + 002) for each atom and defining Q = U2 — u>\. 

There are different physical situations that can be investigated and they are characterized by 
the ratio of the distances between the atoms (r^), which we assume to take fixed positions, and 
the correlation length L c of the reservoir, i.e., jf- << 1 and jf- >> 1. First, let us consider 
that the minimal distance between the atoms is quite large if we compare it to the correlation 
length of the reservoir. All the terms — for all i,j are nearly vanishing and therefore each atom 

interacts with its own reservoir. Since the phase transition from the fluorescent to the superradiant 
emission is a cooperative process involving a collective mode of all the atoms we concentrate our 
investigations in the limit where all the terms r -f- almost vanish. We consider the situation where 
the linear dimension of the total atomic system is small compared to the correlation length L c of 
the reservoir, therefore the N atoms interact collectively with the reservoir. We can proceed by 
assuming that for a fixed /c-mode of the reservoir all the coupling constants go, i — 1, 2, ..., N of 
Eq. (2) are equal. The interaction between the reservoir and the two-level atoms is given by the 
collective pseudo-spin operators J z = Y,f=i a \jy A simplified model is achieved if we assume that 
the two-level atoms and the bosonic modes are in the interior of a high-Q cavity. Assuming that 
the frequency uq of one of the cavity modes is near-resonant with the energy gap Q of the two-level 
atoms, such as situation generates the following physical model: the two-level atoms effectively 
interact only with that mode, and all the other bosonic modes do not couple with the two-level 
atoms. Under these circumstances the model is reduced to a single mode of the bosonic field with 
the creation and annihilation operators tf and b, respectively, interacting with an ensemble of 
atoms. 

Let us define the Fermi raising and lowering operators a\, a i: f3\ and /%, that satisfy the anti- 
commutator relations c^a] + ajaj = 8^ and /3j/3J + PJPi = 5ij. We can also define the following 
bilinear combination of Fermi operators, ajctj — PJPi, a\f3i and finally ftloii. Note that af^, at^ and 



a,-s obey the same commutation relations as the above bilinear combination of Fermi operators. 
Therefore, we can change the pseudo-spin operators of the spin-boson models by using the bilinear 
combination of Grassmann Fermi fields 

°f*) — ► ( a Ui ~ Pi Pi) , (3) 

cr+ — ► a\Pi , (4) 

and finally 



a d) 



P\oa . (5) 



Using the above results, the Euclidean fermionic action for the model that we have defined above, 
can be written as 



Q f A (u*( ^ T K^( *( \®°^Ll^R*( M^ (% Tl < \ 



(6) 



where Hp is given by 



Q N 

H F = u b*(r)b(r) + ^ £(<(tH(t) - A*(r)A(r) 

N 

+4jE{{<*XT)<Xi(T)-ft(TMTJ) (b(r)+b*(r))) . (7) 

Usually, in spin-boson models, the coupling goes with -j=. Since this model has no crossed terms 

of ai and Pi, we are using a -^ dependence in the coupling constant. We interpret such change only 
as a coupling constant renormalization. We can also justify the change adopted using the Holstein- 
Primakoff mapping, where we obtain a divergent Hamiltonian if we use the original coupling -4=. 

With the jj dependence in the coupling constant the Hamiltonian of the system is well behaved. 
Let us define the formal quotient of two functional integrals, i.e., the partition function of the 
interacting model and the partition function of the free model. Therefore, we are interested in to 
calculate the following quantity 

Z_ J[drj\eS 

where S = S(b, b*, a, a\ P, fi^) is the Euclidean action given by Eq. (6), So = 5*0(6, b*, a, afl , P, P^) 
is the free Euclidean action for the free single bosonic mode and the free two-level atoms and 
finally [drj\ is the path integral measure. In Eq. (8) we have functional integrals with respect to 
the complex functions b*(r) and b(r) and Grassmann Fermi fields ot*(r), a^r), P*(t) and Pi{r). 
Since we are using thermal equilibrium boundary conditions in the imaginary time formalism [51] 
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[52] [53] [54], the integration variables in Eq. (8) obey periodic boundary conditions for the Bose 
field, i.e., b((3) = 6(0), and anti-periodic boundary conditions for the Grassmann Fermi fields i.e., 
Oi(P) = -c*i(0) and Pi(0) = -A(0). 

In order to obtain the effective action of the bosonic mode we must integrate over the Grass- 
mann Fermi fields. Therefore, let us define the free action of the bosonic field by 

S (b) = jy T (b*(r)^--u; b*(T)b(T^ . (9) 

The total action of the whole system can now be separated into this free action and a Gaussian 
fermionic part. This is written in the form 

S = S (b)+ drJ2p\(r)M(b*,b) Pi (r), (10) 

where Pi{r) is a column matrix given in terms of Grassmann Fermi fields 

Pi{r) 



Pi( T ) = 

p\(r) = ( i;m n%r) ) (11) 






and the matrix M(b*, b) is given by 

* (m «r))=(* + ' + * WT > +i, < T » n ° |. (12) 

^ o a r - § _ i(6 (T) + (c(t)) ' 

These complex functions and Grassmann Fermi fields b(r), «i(r) and /3j(r) can be represented in 
terms of a Fourier expansion. Therefore, we have 

b(r)=f3- 1/2 J2 b (u)e luJT (13) 

and 

Pt (r)=p- 1 / 2 J2p l (p)e ipT . (14) 

p 

Since the complex function b(r) obeys periodic boundary conditions, and the Grassmann Fermi 
fields OLi{r) and A(r) obey anti-periodic boundary conditions, we have that u> = ^j?- and p = 

t , where they are the bosonic and fermionic Matsubara frequencies respectively. Substituting 



the Fourier expansions given by the Eq. (13) and Eq. (14) in the matrix M pq (b*, b) given by Eq. 
(12) we get 

*,.<O)=( ( * + *>«" + ,. ° \ (15) 



(ip - f)6 pq - Q 



where 



Q = gN- l f3- l/2 (b*{q-p) + b{p-q)\ (16) 

To simplify our calculations, let us change variables in the following way: 

\ 1/2 



7T 



&(«)-> — 6(w) (17) 

\ (CJ - 2Cj) J 

and 

/ \ 1/2 

6»- ? — &». (18) 

Note that Eq. (18) is not the complex conjugate of Eq. (17). It is easy to see that after these 
changes of variables the denominator of the Eq. (8) turns out to be equal to unity 

/[diK&)]e3q>(-7r £&>)&(«)) = 1. (19) 

We can express the ratio J^- by the integral 

|| = /[diri&)]exp(s e// (&)), (20) 

where S e ff(b) is the effective action of the bosonic mode which is given by 

S eff = -ix ]T6*(cj)6(cj) + N\ndet(I + A). (21) 

The matrix A in the determinant of the above equation is given by 

(I + A } J' + E ° ) ,22) 

V I+D J 

where the components of the matrices E and D are given by 



and 



In order to perform the functional integral given by Eq. (20) we must to find a manageable 
expression for det (I + A), so we use the following identity 



z 
Z 



det (I + A) = det (I + E+ D + ED) ~ e *E+trD-\E*-\D* ^5) 

where the approximation in the exponent on the right-hand side of Eq. (25) is performed up to 
second order in the E and D matrices. Therefore, upon using the expressions of E and D given 
in Eq. (23) and Eq. (24), respectively, into the approximation given in Eq. (25), we can find an 
expression for the ratio -^- defined in Eq. (20) of the following form 

[drj(b)\ exp (-7T Y, b *(u) b(u) + a (6(0) + b* {0)] 
- I E c(uO (&M K-co) + 6» b*(-u)) - 1 £ d(u) (b(co) &») V (26) 

where the terms ao, c{uj) and d{ui) are, respectively, defined in the form 



M p(<*o - *») r (p 2 + f ) ((? - -) 2 + f ) • ( j 



It is possible to show that the contributions c(cj) and d(u) cancel out for u ^ and c(0) and d(0) 
are different from zero, but their values do not matter in the discussion next because in the limit 
N — > oo they do not contribute to the partition function. Therefore in the thermodynamic limit 
N — ► oo we have that -§- can be written as 

■f- = / I] ^)^(w)e" ,E ^ l ' MW /rf6(0)^*(0)e-^*( ) 6 W +ao ( 6 (°) +fe *( )) . (30) 



The integrations for uj 7^ are not dependent on thermodynamical parameters. Thus it can be 
considered as a normalization constant, i.e., 

-J- = C I db(0)db*(0)e^ b *^ b ^ +a »( b W +b *^ . (31) 

The partition function is therefore an entire function. From the In Z quantity, which is given by 

lnZ = lnZ + lnC , + ^tanh 2 ^V (32) 

ujq V 4 / 

we can find the free energy of the model. Also the mean energy E and canonical entropy S can 
be computed. We get 

E = E °-^ tanh {-) eosh 2 (M) ' ^ 



qi&n tanh(^) 

S = S + \nC - 9 -£- -^( (31) 

2 ^0 cosh 2 1 iJ f 



and 



where So is the canonical entropy for the free case. Using the third law of thermodynamics we get 

C = l. 

We can see that the contribution of the last term in the right-hand side of Eq. (32) is coming 
from the zero mode. Owing to the absence of the coupling via at--, and aV--, to the bosonic mode, the 
zero mode generates a negative term in the expressions for the energy and the entropy. It is known 
in the literature that the zero modes can lead to problematic results [55] . In the models discussed 
by Dowker, the entropy becomes negative and also a temperature-dependent pole appears in the 
free energy [56]. At this point, as discussed before, there are two different paths that we can 
follow concerning the zero-mode problem. One is to disregard the zero mode, thus solving the 
problem of the sign of the canonical entropy. We would like to stress that we can not identify 
the zero-mode as a Nambu-Goldstone mode [57] [58] since there is no spontaneous breaking of 
the continuous symmetry in the model. The second path is to assume that g 2 has some bound. 
For small values of g 2 the model has positive entropy for finite temperature and it goes to zero 
at /3 — > 00, satisfying the third law of thermodynamics. Both procedures, the ad hoc first one of 
throwing away the zero-mode, or assuming some bound in the g 2 value, solve the problem of the 
negative entropy, making the model consistent from the thermodynamical point of view. Adopting 
any of these procedures, the partition function is analytic for all temperatures. 



3 The functional integral for the generalized fermionic 
Dicke model 

The generalized Dicke model, where an ensemble of identical iV atoms interacts linearly with 
one mode of a bosonic field, is defined by the Hamiltonian 

JV q N 

H = h ® E ■w<?lj)+wotfb®lB+ -y= E( 6 + 6t ) <8) ( t7 5) + %))- ^ 

j=l z ViV i=1 

In the above equation g is the coupling constant between the atom and the single mode of the 
bosonic field. The b and 6' are the boson annihilation and creation operators of mode excitations 
that satisfy the usual commutation relation rules. 

The aim of this section is to prove that a model with an interaction Hamiltonian generating 
only virtual processes presents a phase transition from normal to superradiant state at some 
temperature with the presence of a condensate and also a quantum phase transition at some 
critical coupling. Similarly to the last section we find that the Euclidean action S is given by Eq. 
(6), where now Hp is the full Hamiltonian for the generalized fermionic Dicke model given by 

Q N / \ 

H F = u b*(r)b(r) + - E (o.-(T)oi(T) - #(r)ft(r)J + 

Ui E k*(T)&(T)6(T) +ai (T)(3*(T)b*(T)) + 



Vn ti 



9i 



N 



+ 77= E («,(r)/3*(r)6(r) + a*(r) &(t) b*(r) ). (36) 



Note we are introducing two coupling constants, g\ and g 2 , for the rotating and counter- rotating 
terms, respectively. As we discussed before, the main reason for this is that we are interested 
in to identify the contribution of the real and virtual processes in the phase transition with the 
formation of the condensate. We are interested in calculating the formal quotient given by Eq. 
(8). Following the same calculation from last section we can arrive at Eq. (20) and Eq. (21), 
where the matrix A is now given by 

In the equation above the quantities B and C are matrices with components given by 

i _i i 

7r\ 2 /. fi\ 2 ( g\b*{q — p) gib (p — q) \ (. Q 



B >« - [W) [ ip + 1) [ ^-ib-r) + ^-.-(p-,) 1 ' k - - ' (M 
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and 



a 



/"/ 



7T 



(3N 



ip 



rr 



9ib(p-q) + g 2 b*(g-p) \ / + ft\ 5 
w -i(p-g) Jw -i{q-p)/ V 2 / 



(39) 



We shall investigate the integral given by Eq. (20) for temperatures that satisfy /3 x > f3 c 1 . First 
of all let us show that this integral converges. Using the estimate 



| det(J + A) | < exp ( Re (trA) + - tr(AA ] ) 
we can show that the ratio -^- obeys the following inequality 



(40) 



Z 
— < 



n 

w>0 



- oo(0) + 2c (0)) (l - o (0) - 2c (0) 
( 1 - o (o;) + 2c (^)j f 1 - o (o;) - 2c (cj) J 



-i 



(41) 



where the cio(uj) and co(cj) are given respectively by 



a {uj) 



g\ + g\ 



E 



i 



i 



and 



Cq(u) 



0(u% + W 2 )V2 p _^ w { f + g 2)l/2 (^ + p 2)l/2' 

1 1 



^ogig2 



E 



(42) 



(43) 



/?K 2 + ^ 2 ) p _v =w (¥ + ? 2 ) 1/2 (¥ + p 2 ) 1/2 ' 

In a similar way like Popov and Fedotov [28] proved, for the case of rotating wave approximation 
we have that < a (uj) + 2c (u) < a (0) + 2c (0) and a (0) + 2c (0) = 0(c<j~ 2 ln a;). Therefore 
if a (0) + 2c (0) < 1, then Eq. (41) guarantees convergence of the expression -^-. The condition 
<2o(0) + 2cq(0) = 1 is the equation for the transition temperature, then we have 



a (0) + 2c (0) 



gi + g2 

f2o;o 



-(¥ 



1. 



The inverse of the critical temperature f3 c is given by 



4 tapir 1 f QU ° 



(44) 



(45) 



Note that there is a quantum phase transition where the coupling constants g\ and g 2 satisfy 
g\ + #2 — (^o ^) 5 • For ji + j2 / (^o ^) 5 the partition function is no more an entire function in 
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the positive half of the complex plane for the temperature P~ x given by Eq. (45). The system 
enters in a superradiant phase. To calculate the asymptotic behavior of the functional integrals 
at temperatures that satisfy /3" 1 > (3" 1 , we can do the following approximation 



A? 



det N (I + A) = det N (I + BC) -> exp (Ntr(BC) . 



(46) 



After some calculations -§- can be written as 

^0 



Z 

~Z~o 



n 

w>0 L 



- o(0) + 2c(0)) (l - o(0) - 2c(0) 
(l - 0(0;)] (l - a(-uj)) - c 2 (uj) 



■ 1/2 



n-i 



+ 



+0(N- 1 ), 
where a{ui) and c(u>) in the above equation are given, respectively by 



a{uj) 



and 



c(u) 



' g\ (Q - iu) _1 + g\ (Q + iu>) 
(w - iu) 

( 9i92& N 



-r 



tanh 



pn 



tanh 



pn 



(47) 



(4£ 



(49) 



V(cu 2 + cu 2 ) 1 / 2 ^ 2 + CJ 2 ) 

Taking the limit (N — * 00) in Eq. (47) we are in the thermodynamic limit. We turn out to the 
discussion concerning the local elementary excitation of the ground state. To find the collective 
excitation energy level spectrum we have to use the equation 



c 2 M 



0(0;) 



1 — of — uj) 



0, 



(50) 



and making the analytic continuation (iu —> E), we obtain the following equation 



+ 



= — 


4 1 4 

91 + 91 
_{ul - E 2 )(Q 2 - E 2 )_ 


tanh 2 f-— ") + 






[ 9\9l ( 1 1 4fi 2 y 


u*-e?) 


+ 


>g - E 2 )\(n -E) 2 (n + E) 2 (tt 2 -E 2 ) 2 J_ 


\gK 


Q - E)- 1 + g 2 {Q + L 
(u - E) 


J)- 1 i gfttl + E)- 1 +g 

(u + L 


2 (Q- E)- 1 - 


-(*? 



(51) 
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Solving the above equation for the case (3 l = f3 c 1 we find the following roots 

E 1 = (52) 

and 

E 2 = ( gi^ + "o) 2 + ^(^-^o) 2 \ V2 (53) 

V (91 + 92) J 

Its low energy state of excitation is a Nambu-Goldstone mode. Now, let us present the critical 
temperature and the energy level spectrum of the collective bosonic excitations of the model with 
the rotating-wave approximation, where g\ 7^ and g 2 = 0. The result obtained by Popov and 
Fedotov is recovered, where the equation 



and 

9l 



o(0) = 1 (54) 



tanhf^— )=1, (55) 



give the inverse of the critical temperature, f3 c . It is given by 

P c = — arctanh ( — 5- ) . (56) 

" V 9\ J 

The order parameter of the transition is the expectation value of the number of excitation associ- 
ated to the bosonic mode per atom, i.e., lirriN ^oc ^j^- 7^ 0. Note that again u , £1 and g\ define 
also a non-zero critical temperature where the partition function is no more an entire function in 
the positive half of the complex /3-plane. We may expect a superradiant phase for the temperature 
flc 1 given by Eq. (56). The energy level spectrum of the collective Bose excitations in this case is 

£i = 0, (57) 

and 

£ 2 = fi + cj • (58) 

In this case, there is also a quantum phase transition, i.e., a zero temperature phase transition 
when gi = (ujq Q)~ . Now we will show that is possible to have a condensate with superradiance in 
a system of N two-level atoms coupled with one mode of a Bose field where only virtual processes 
contribute. In the pure counter- rotating wave case, i.e., g\ — and gi 7^ 0, the inverse of the 
critical temperature, f3 c is given by 

P c = — arctanh ( — 5- ) , (59) 

Q V gi J 
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and the spectrum of the collective Bose excitations given by 

#i=0, (60) 

and 

E2 = \n-u \. (61) 

A comment is in order concerning the spectrum of the Bose excitations. In both of the cases: 
working with the pure counter-rotating or the rotating-wave terms, there is a phase transition. In 
the case of the rotating-wave approximation gi ^ and gi = 0, there is a Nambu-Goldstone mode 
(E = 0). In the pure counter-rotating case gi = and g 2 ^ also there is a Nambu-Goldstone 
(gapless) mode. Thence we show that it is possible to have a condensate with superradiance in a 
system of iV two-level atoms coupled with one mode of a Bose field where only virtual processes 
contribute. Since the energy E 2 ~ 0, local elementary excitations of the ground state with low 
energy can easily be created causing a significant fluctuation effect. Unfortunately we are not able 
to evaluate these effect in the systems. 

An important question is the way of practical realization of the second model, i.e., the gener- 
alized Dicke model in the laboratory. As was stressed by Dimer et al [49] it remains a challenge 
to provide a physical system where the counter- rotating terms are dominant. Experimental ob- 
servation of the superradiant phase in a situation where is possible to control the importance of 
the counter-rotating terms in the generalized Dicke model (g\ ~ 0, g% ^ 0) could improve our 
understanding of this phenomenon. 



4 The functional integral for the fermionic Dicke model 
with intensity-dependent coupling 

A model where the behavior is quite interesting from the physical point of view is the one where 
the coupling between N identical two-level atoms and one mode of a bosonic field in a lossless 
cavity is intensity-dependent. The generalization for this model introducing the counter-rotating 
terms is straightforward. Therefore, in this paper we discuss the model with rotating and counter- 
rotating terms, where the intensity-dependent contribution appears only in the real processes. The 
necessity for disregard the intensity dependent coupling contribution in the part of the interaction 
Hamiltonian that generates the virtual processes is easy to justify. It is counterintuitive that the 
virtual processes are amplified when the number of excitations in the bosonic sector is enhanced. 
We will show that at low temperatures the contribution from the counter-rotating terms dominate 
over the rotating ones and the system presents a quantum phase transition. 

After this discussion, the total Hamiltonian of the atoms and the bosonic mode takes the form 

N Q 

H = I s ® E 17 a (j) + ^bH ® I B + 



^ 2 
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N 

Y, (b (b j b)% <g> a^ + b ] (b j b)? <g> (Ty } + & f <g> (tJ } + & <g> cr^) . (62) 

At this point, let us make a parallel between the model given by Eq. (62) assuming the rotating- 
wave approximation, and the Jaynes-Cummings model [42], where one two-level atom is coupled 
to a single mode quantized electromagnetic field. It is possible to generalize this model, by using 
the rotating- wave approximation, known as the two-photon Jaynes-Cummings model [59] [60]. 
The introduction of the counter-rotating terms is straightforward. It has been discussed in the 
literature the possibility of controlling the relative importance of them in laboratory using an ion 
trap [50]. 

Going back to the model we have also non-Gaussian functional integrals to solve. The main 
point we wish to demonstrate in this section is that the model defined by the equation bellow 
at low temperatures is exactly soluble and the partition function can be presented in a closed 
form. Changing the atomic pseudo-spin operators by a linear combination of Grassmann Fermi 
field yields the fermionic generalized Dicke model with intensity-dependent coupling defined by 
the Hamiltonian 



V ^ 

2 



H = Eo*0-)+<"b*t6- 



i=i 



N N 



-%J2( b ( &t &)*«& + bj ( &t &)*/&*) + -7= E (*#«* + &t «!a) , (63) 

ViV j=i ViV j=1 

where again we are introducing different coupling between the single mode bosonic field and the 
reservoir, g\ and #2 for rotating and counter-rotating terms, respectively. The difference of the 
Hamiltonian for this model, given by Eq. (63), when contrasted with the Hamiltonian of the 

generalized Dicke model, Eq. (36), is the term (tf 6) "a in the part of the interaction Hamiltonian 
which generates the virtual processes. The fermion integration on the generalized Dicke model 
is repeated for this case, so we can take expressions given by Eq. (20) and Eq. (21) where the 
matrix A is given by 

In the equation above the quantities B and C are matrices with components given now by 

Dpq — 





,|,,i — 1 (r(0)K0))* , b ' iq : P) +«, , b{p - q) I «w> 
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and 



N/3J V 2 7 V 2 . 

*Ulf| (6-(Q)KQ))' . t(p - g) + », , b ' (q ' P) I. (66) 



For this case, in order to find an approximate expression for the partition function in the ther- 
modynamic limit (N — ► oo), we use the approximation det(J + A) = det(J + BC) ~ exp(trBC). 
Similar procedure was used for the generalized Dicke model case. Then we can arrive at the 
following approximate expression 



|- = J[d V (b)]e X p(-nJ2{±-c(uj))b*(uj)b(u) 

+ nY, a M b *(°) b (°) b *( u ) 6 M + *" E d ( u ) ( b *(°) b (° 

xfb(to)b(-uj) + b*(u)b*(-u) H , 
where a (a;) and c(o;) of above equation are given, respectively, by 



(67) 



a(a;) =- — tanh — — — —— : — -, (6l 



c(u) = g\ tanh( — - ) — , (69) 



/3fi\ 1 

(Q + iuo) (ujq — iu) ' 
and 

d(w) = ai0 2 [7^- ) tanhf-— ) r . (70) 

1 ; yiy2 \[3u ) \ 4 )(n + iuj)(u% + cu 2 )5 l ; 

In order to perform the functional integral given by Eq. (67) we have separated the zero and 
non-zero modes, so we can write 

* = /^(0)^(0)e< a(0)( " ( ° ))V(0) + (c(0) - 1)6(0) " (0) ) 

Zq J 

i E^( 1 - c ( u )) l, W( u ) 



I X\db(uj)db*(u)e~^^ 



xe » E^o a ^ b *(°) fe (°) 6 *H 6M + d(w) (b*(0) 6(0)) 5 (b(w) 6( -«) + 6*(w) b*(- w) ) /y-Q 
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It is quite difficult to compute in a close form the quantity -J^- since the integrals that appear in 

Eq. (71) are not Gaussian. Nevertheless, at low temperature (/3 — > oo), we can make it evident 
that the contribution coming from the counter-rotating terms dominates over the rotating ones. 
We can take this limit in the Eq. (68), Eq. (69) and Eq. (70) so we get 

lim — = / ' n^)d&V)e~ wS « (1 ~ c(w))6 * (w)6(w) - (72) 

Performing this last integral we have 

lim ^- = U(l-c(u))~\ (73) 

where the term c(u) in the limit (3 — ► oo is 

lim c{u) = f — . (74) 

/3^oo (i I -f too) (uj — lUU) 



i 



From this last equation we can see that the ratio Z/Zq is non-analytic for gi — (o>of£)2, thus 
corresponding to a quantum phase transition. Therefore, in the model at some range of the 
physical parameters, the counter-rotating terms are dominant, leading to a phase transition. Since 
we are not able to evaluate the ratio ^- given by Eq. (71) for all temperatures, it is not possible to 
describe the phase diagram of the system near the quantum phase transition. The possibility that 
appears a line of second-order phase transition for /3 _1 > terminating at the quantum critical 
point must be investigated. 



5 Conclusions 

In the present paper we are considering the issue of the formation of the condensate with a su- 
perradiant phase transition in spin-boson models with a single bosonic mode with quite particular 
couplings between the single bosonic mode and a environment of N two-level atoms. We have 
assumed that the whole system is in thermal equilibrium with a reservoir at temperature f3~ l . The 
interaction Hamiltonian of the models also generate virtual processes. We first investigate a model 
based in one discussed by Chang and Chakravarty, Legget and others. Second, with the couplings 
gi and g 2 for rotating and counter-rotating terms, respectively, we has defined the generalized 
Dicke model. Finally we study the generalized Dicke model with an intensity-dependent coupling. 
In this last model it is necessary to use perturbation theory to investigate the thermodynamic of 
the model. Nevertheless in the low temperature limit the model is exactly soluble. 

Studying the case where in identical two- level atoms act as a thermal reservoir (N — » oo), 
we investigated the thermodynamics of the three models using the path integral approach with 
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functional integration method. For generic g 2 , we found that the first model is unrealistic from 
the thermodynamical point of view, since the entropy becomes negative at some temperature. 
As we discussed, there are two different solutions for this problem. The first one is an ad hoc 
procedure, disregarding the zero-mode, that can make the model admissible. Another possibility 
is to assume that g 2 is a bounded quantity. After any of these procedures, it is shown that the 
partition function is analytic for all temperatures, consequently there is no phase transition in 
the model. In the second model the situation is more interesting. We studied the nonanalytic 
behavior of thermodynamical quantities of the generalized Dicke model by evaluating the critical 
transition temperature and presenting the spectrum of the collective bosonic excitations, for the 
case g\ 7^ and gi = 0, g\ = and g-i 7^ and also in the general case. Our result show that it is 
possible to have a condensate with superradiance in a system of N two-level atoms coupled with 
one mode of a bosonic field where only virtual processes contribute. It is important to realize that 
the energy of the mode, which is not the Nambu-Goldstone mode in Eq. (58), is always larger 
than the one in Eq. (61), i.e., in the system where the condensate appears due to the virtual 
processes. Therefore both processes, real and virtual ones give different contributions to generate 
the condensate. In the generalized Dicke model with the intensity-dependent coupling we found a 
superradiant phase transition at zero temperature, i.e., a quantum phase transition. 

One interesting aspect of superradiance is the fact that initial quantum fluctuations that are 
able to cause spontaneous emission in a few atoms are amplified leading to fluctuations on a 
macroscopic scale. We would like to point out that another mechanism to enhance the importance 
of virtual processes, in quantum field theory in the presence of macroscopic boundaries, was 
proposed by Ford [61]. In a local calculation, Ford presented a mechanism based on the fact that 
the contribution of various part of the frequency spectrum of the Casimir effect is an oscillatory 
function. The contributions of different ranges of frequency almost, but not quite, cancel out 
one another, and there is the possibility of enhancing the magnitude of the effect by altering 
the reflectivity of the boundary in selected frequency ranges. A quite different mechanism for 
amplification of vacuum fluctuations was proposed by Ford and Svaiter [62] [63] using parabolic 
mirrors to produce large vacuum fluctuations near mirror's focus. These authors studied the 
renormalized vacuum fluctuations associated with a scalar and electromagnetic field near the 
focus of a parabolic mirror. Using the geometric optics approximation they found that the mirror 
geometry can produce large vacuum fluctuations near the focus, similar to what happens in the 
classical focusing effect by the parabolic mirror geometry. 

With the development of quantum information and its application in computation and com- 
munication the entangled states have been attracting enormous interest, since several quantum 
protocols can be realized exclusively with the help of entangled states. In the spin-boson model 
with the dipole-dipole interaction included, one can address the issue of how does the dipole-dipole 
interaction change the critical temperature, where again the system exhibits a phase transition 
from the fluorescent to the superradiant phase. An important subject for the near future is to 
study the degree of entanglement between the N two-level atoms and the bosonic mode near the 
quantum phase transition in the generalized Dicke model [64] [65] [66]. This subject is under 
investigation by the authors. 
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